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Abstract 

In a recent FTC by Tremblay et al (2009 J. Phys. A: Math. Theor. 42 205206), 
it has been conjectured that for any integer value of k, some novel exactly solvable 
and integrable quantum Hamiltonian Hk on a plane is superintegrable and that the 
additional integral of motion is a 2/cth-order differential operator Y^fc- Here we demon- 
strate the conjecture for the infinite family of Hamiltonians Hk with odd k > 3, whose 
first member corresponds to the three-body Calogero-Marchioro-Wolfes model after 
elimination of the centre-of-mass motion. Our approach is based on the construction 
of some -D2fc- ex te n ded and invariant Hamiltonian Hk, which can be interpreted as a 
modified boson oscillator Hamiltonian. The latter is then shown to possess a Dik- 
invariant integral of motion 3^2fc) from which l^fc can be obtained by projection in the 
D2k identity representation space. 
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1 Introduction 



In a recent work, Tremblay, Turbiner and Winternitz (TTW) have introduced an infinite 
family of exactly solvable quantum Hamiltonians on a plane 

H k = -d 2 T - -d r - \d 2 + co 2 r 2 + ^(asec 2 kip + (3 esc 2 kip), (1.1) 

defined in terms of polar coordinates (r, (p) and of an arbitrary real number k pQ. In (11. II) . 
the configuration space is given by the sector 0<r<oo,0<<y9< n/(2k), and u, a, (3 are 
three parameters such that u > 0, a = a(a — 1) > — l/(4/c 2 ), = b(b— 1) > — l/(4k 2 ). The 
Hamiltonian H k has been shown to be integrable with a second-order integral of motion 

X k = -&l + k 2 {a sec 2 k<p + P esc 2 k<p), (1.2) 

responsible for the separation of variables of the corresponding Schrodinger equation in 
polar coordinates. 

In the same work, it has also been conjectured that for all integer values of k, the 
Hamiltonian (11. ip is superintegrable and that the additional functionally-independent in- 
tegral of motion is a 2fcth-order differential operator Y^k- This conjecture has been proved 
by a direct calculation of Y 2 k for k — 1, 2, 3 and 4, as well as for k — 1 to 6 and 8 in 
the special case where either a or vanishes. Actually, for k — 1, 2 and 3, the super- 
integrability of H k was already known since the Hamiltonian then reduces to that of the 
Smorodinsky- Winternitz system (k = 1) [21 E], of the rational BC2 model (k — 2) [1] or of 
the Calogero-Marchioro-Wolfes (CMW) model (k = 3) 0(6]. 

More recently, the credibility of the conjecture has been reinforced by showing that for 
the corresponding classical systems all bounded trajectories are closed and that the motion 
is periodic for all integer and rational values of k [7] . Such classical systems have then been 
proved to be superintegrable [S] and generalizable to higher dimensions [9]. 

Since it is clear that the proof of the conjecture for k > 4 requires a different ap- 
proach from that based on a brute-force calculation of Y 2 k, it is worth investigating the 
use of some differential-difference operators or covariant derivatives, known in the math- 
ematical literature as Dunkl operators [10]. Such a formalism, based on the introduction 
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of exchange operators, has indeed been successful in proving the superintegrability of the 
iV-body Calogero model [Til [12], which for N = 3 is related to if 3 with (3 = 0, as well as 
that of the CMW model [T3] , for which some results in polar coordinates are also available 



In a first step along these lines, we have generalized the exchange operators known for 
the CMW model to the whole set of Hamiltonians H k with integer k [12]. For this purpose, 
we have realized the elements of the dihedral group D 2k as operators on the plane (r, <p) and 
employed them to define some differential- difference operators D r and D v , where the latter 
assumes a different form according to whether k is odd or even. We have then constructed 
/^-extended and invariant Hamiltonians from which the starting Hamiltonians H k can 
be retrieved by projection in the D 2k identity representation space. Furthermore, we have 
shown that X k can be obtained by the same procedure from the operator —D^, commuting 
with % k . 

The purpose of the present communication is to build on these advances in order to 
provide a rigorous proof of the above-mentioned conjecture for any odd k value. More 
specifically, we plan to demonstrate the existence of a second integral of motion Y 2k of H k , 
obtainable by projection from a corresponding operator y 2k for the extended Hamiltonian 
Hk, whenever k is odd and greater than or equal to three. 

2 Exchange operator formalism for Hk with odd k 

In [T5], it has been shown that for any integer k value, the 4k elements of the dihedral 
group D 2 k can be realized as W and 1Z 1 T, i = 0, 1, . . . , 2k — 1, where 1Z = exp (kird^) is 
the rotation operator through angle ir/k in the plane (r, cp) while X = exp(iircpd v )2 changes 
(p into —(p. Such operators indeed satisfy the defining relations of (a unitary representation 
of) D 2k , 

U 2k = j2 = 1; xr = n™- 1 !, -r) = n 2k -\ X f = X. (2.1) 



1 In the definition of 1 and only there, i stands for the square root of — 1. 
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The partial derivatives d r and d^ can then be extended into differential-difference oper- 
ators D r and D 9 , defined for any odd k by 

'k-i 



D r = d r - -{aK + b) ( n2i j X > 



,i=0 



k-l 



k-l 



Dtp = d„ + a^tan ( cp + ^ ) U k+2i X - 6^cot [<P+j) T^'l. 



(2.2) 



i=0 v ' i=0 

where a and b are related to the potential strengths a and /3 in H^- From (11. ip . it follows 
that D r and .D^ fulfil the following Hermiticity, exchange and commutation relations 



D\ = -D r -- 

r 



l + 2(aK + b) (^K 2 ^! 



Dl 



-D 



UD r = D r TZ, XD r = D r I, UD^ = D V K } 
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'k-l 



XD U 



(2.3) 



[D r ,D v ] = —(aK + b) ^2n 2i )XD v , 



(2.4) 



i=0 



respectively. It is worth noting that the first relation in (I2.3P as well as equation (12. 4p 
significantly differ from the corresponding results for d r and d v , namely dl = —d r — - and 
[d r , d v ] = 0. 

From D r and D v , one can build a /^-extended Hamiltonian 

l + 2(aU + b) (j^ 71 } 1 



U k = -Dl - - 



,i=0 



D r --Dl + uj 2 r\ 



(2.5) 



which can also be rewritten as 



Hi 



-dl - -d r - 



k-l 



D% ~k{a 2 + b 2 + 2abK) ^ TZ 2i 



i=0 



i 2 2 

+ u r 



with 



k-l 



k-l 



i=0 



D l = d l-Y, sec2 [v + i) < a - nk+2l% ) - E csc2 ( v + t ) h{h ~ n2lJ) 

k-l 

+ k{a 2 + b 2 + 2abn)^K 2i . 



i=0 



Such a Hk is clearly invariant under D^k and, due to some simple trigonometric identites, 
gives rise to after projection in the D 2 k identity representation space, i.e., after replacing 
1Z and I by 1. 



Furthermore, X k = —D 2 is a D2k invariant, commuting with Hk, and whose projection 
in the same space is given by X k — k 2 (a + b) 2 , with X k defined in equation (II. 2p . 



3 Modified boson operators for odd k > 3 

To go further, it is worth observing that, on one hand, the starting Hamiltonian H k has 
an oscillator spectrum [I] and, on the other hand, the extended Hamiltonian defined 
in (I2.5p . has (apart from some additional exchange operators) a form reminiscent of that 
of an oscillator Hamiltonian in polar coordinates, h = —d 2 — ^d r — ^.d 2 ^ + u 2 r 2 . Since the 
latter can also be written as h = 2u(h a + a&), h a = |{at, a}, = \{b\ b}, in terms of two 
sets of boson creation and annihilation operators a', a and fet, b, it seems appropriate to 
introduce for %k two sets of modified boson operators 



and 



.4 
A* 

B 
fit 



2u 
1 



1 



2u 
1 



cos (p(cur + D r ) sin ipD^ 



cos (piujr — D r ) H — sin <pDu 
r 



sin <p(ujr + D r ) H — cos tpD^ 



sin if(ur — D r ) cos <pD^ 



(3.1) 



(3.2) 



by simply substituting D r and for d r and d v in the definition of a, at, b, &t in polar 
coordinates. 

At this stage, it should be noted that the fact that and B^ are the Hermitian 
conjugates of A and B, respectively, is not so obvious as for at and 6t with respect to a 
and 6. To prove such a property, one has actually to use equations fl2.ll) . (12. 3p and the 
definition of 1Z and X, as well as the commutation relations 



[D r ,r] = l, 
k-i 



[D v ,r) = 0, 



[D r , cos <p] 



i=0 



i=0 



sin ( cp + — j sin — K 2l X 



(3.3) 



(3.4) 



[D r , sin ip] 



k-l 



—2a cos I if + 



i=0 



VK 



k 



k-l 



i=0 



+ 26 2^ sm I + — J cos —Tl 2% X 



(3.5) 



k-l 



\. — \ / 27T \ * 7r /ni--i-2»'-T- 

[TJ^, cos ip\ = — simp — 2a > sin I ip + — 1 cos —7c ^ I 



i=0 



fc-i 



t=0 



+ 26 2^ cos I ip + — 1 sm —T&T 



(3.6) 



fc-i 



[D v , sin ip] = cos ip + 2a ^ sin ( ip + — J sin —TZ k+2i X 



i=0 



k-l 



8=0 



26 2^ cos ( + — ) cos —7c 2l l 



(3.7) 



which can be easily established from (12. ip and (12. 2p . 

It turns out, however, that the two sets of modified boson operators (13. ip and (13. 2p 
do not have simple transformation properties under D2k whenever k > 1. To cope with 
this drawback, it is advantageous to replace A, A* by k dependent pairs of modified boson 
operators A iy A\, i = 0, 1, . . . , k — 1, defined by 



Ai 
At 



2oo 
1 



ITT 



cos ( (p + — ) (ur + D r ) — — sin ( <p + 
cos [<p + l — ] (ur — D r ) H — sin ( (p + 



D„ 



D u 



(3. 



and such that Aq = A and A\ = A^ . In the same way, B, may be replaced by Bi, 
Bj, i — 0, 1, . . . , k — 1, obtained from (13.21) by substituting <p> + ?f for ip (hence B = B, 
Bj = Bt). For most purposes, though, it will be enough to only consider Ai, A], since Bi, 
B\ can be expressed in terms of them as^| 



fc-i 



sm 



(3 ~ J)7T 



fc-1 



sm • 



i=o 



(j ~ J> 



(3.9) 



2 Up to equation (|3.8p . all mentioned results are valid for any odd k value, including k = 1. This is not 
the case for equation (|3.9[) . nor for some formulas to be given in the sequel, which require k > 3. From now 
one, we therefore restrict ourselves to such values. 
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Under Z?2fc, the operators Aj, i — 0, 1, . . . , k — 1, transform among themselves (and 
similarly for their Hermitian conjugates). We indeed obtain 

ftA^ 1 = A i+1 , z' = 0,l,...,/c-2, TZAk-iW 1 = -A , 

(3.10) 

XAoX-^Ao, lAfl- 1 = -A k _ i} i = 1,2,..., A; -1. 

These relations can actually be simplified by extending the definition of Aj (and A\) to all 
values z 6 Z, 

A Afe+l = (-1) A A,, z = 0,l,...,A;-l, AGZ. (3.11) 

We then get 

KA{RT X = A i+1 , lAil- 1 = -Ak-i, % e Z. (3.12) 

Among the operators Aj, i = 0, 1, . . . , k — 1, there are only two linearly independent 
ones (e.g., Ao = A and i?o = B). The Aj's therefore satisfy k — 2 independent linear 
relations. We quote here some of them, 
fc-i fc-i 

A: nns 

k 2 



^(-l) l A, = 0, ^A iC os^ = |a , (3.13) 



i=0 i=0 

Z7T I 27T 1 

A.COS— = -(A + A 2i ), fi.sin— = -(A - A 2i ), z 6 Z, (3.14) 

which are easily derived through elementary trigonometry and will subsequently play an 
important role. 

The commutation relations of Aj, Aj, i = 0, 1, . . . , k — 1, can be determined from 
equations (12 .4p and ( 13. 3p . as well as from the relations obtained from H3.4[) - f !3.Tj) by acting 
with on both sides. After a rather lengthy, but straightforward calculation, we get the 
simple results 

[Aj, Aj] = [Aj, A]] = 0, (3.15) 
[A,, Aj] = [A,, Aj] = cos + 2a £ cos cos td)l n ^ 

rv rv rv 

, (l — i)n (I — ?)7T 9 , 

+ 26 > sin v - 7 sin ^ — -^-K 2l l, 
^ k k 

i 

for z', j = 0, 1, . . . , k — 1. Note that from now on, all summations over I run from to k — 1. 
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4 Modified boson oscillator Hamiltonians and super- 
integrability of for odd k > 3 

The modified boson operators A|, i = 0, 1, k — 1, can be used to define some 
modified oscillator Hamiltonians as 



i = 0, I. 



From (jjnOJ), it follows that 

IZHiJZ 1 = ii/j+i, 



0,1,..., A; -2, KHk-iK- 1 = H , 



XHqX 1 = i? , XifjX 1 = H k _i, i — 1, 2, . . . , k — 1, 



(4.1) 



so that the ifj's transform among themselves under Z^fc. 

The explicit expression of in terms of r, (p, D r and -D v can be easily found by first 
considering % — and then acting with 7£ l on the result. It is given by 



2ooH i 



COS" 

1 



(ITT \ ! / %7l \ ( 17C \ 

up + — \D 2 r + - sinf p + — J cosf <p + — j (D r D v + D v D r ] 



^ Sil1 [ (p+ Jl D f 



sm" + — j + 2a ^2 



cos 2 (_L^}l n k + 2i x + 2b y sin 2 il^3l n ^x 



I 



k: 



D r 



+ 



27T \ / ITT 

— 2 sin [ y? + -T- 1 cos I <£> + — 



2"£ 



sm 



(Z — i)7r (l-i) 



7T 



COS ■ 



A' 



, (1 — 1)71 (I — i)7T 9 , 

+ 26 > sin v , ; cos V , ; ^ 2 'X 



+ tu 2 r 2 cos 2 I y + 



I7i 



Elementary trigonometry now shows that the /^-extended Hamiltonian Tik of equation 
f)2.5p is directly connected with the H^s through the relation 

fc-i 



k 



2uj 2J #i = ^^fc- 



(4.2) 



i=0 



Hence, it may be regarded as a modified boson oscillator Hamiltonian. 

Next we may ask ourselves what are the commutation relations of the H^s. As explained 
below, for our purpose it is actually enough to determine the commutator of H with 
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z = 1,2, 



k — 1. On using the convention (13. lip for simplicity's sake, as well as equations 



(I3.12p . ( I3.15P and (I3.16p . we get after some calculations 



[H , H] = (A Ai - A\A ) cos j + 2a ^ 



-{A\A 



, t , . Itt (I - i)-K 
2k + 2l-i + A\A k+ 2l) cos — cos 



+ -(AlA -A{ +2l A k+2l ) 



cos 



(Z-i)tt 



k 



+ l(Al k+2l ^A 2k+2l -i - A\Ai) cos 2 



Zvr 



^ fc+2 'z 



(Z - z)tt 

ife 



+ 26^ -(4^fe+2«-i + A A 2i) sin y sin ^ + ^{A\A - A ] 2l A 2l ) sin 

+ \{A\ +2l _ % A k+2 ^ - A\ A,) sin 2 £j W»I - | sin ^ ^ sin " + 62) ^ 2 '- 

(4.3) 



It may be observed that this relation is also valid for % = 0, leading to [H , H ] = as it 
should be. 

On summing both sides of equation (14. 3p over % from to k — 1 and employing equations 
(13.131) . (13.141) . as well as some relations easily obtained from them, we arrive at the result 
[H ,7i k ] = 0, where we have taken equation (14.21) into account. It only remains to apply 
1Z 1 and to use the invariance of % k under D 2k , together with equation (14. ip . to infer that 



[Hi, n k ] = o, 



i = 0, 1, . . . , k — 1, 



showing that Hi, i = 0, 1, k — 1, are integrals of motion of H k . 

From them, we can form two D 2k invariants, namely their sum proportional to H k and 
their symmetrized product 



y 2k = (2u) k 2J Hp(o)H p (i) ■ ■ - H p 



(k-i), 



(4.4) 



where the summation runs over all k\ permutations of 0, 1, . . . , k — 1. Projection in the D 2k 
identity representation space then leads to H k , on one hand, and to an integral of motion 
Y 2k of the latter, on the other hand. Since Y 2k is clearly a differential operator of order 2k, 
the conjecture of pQ will be proved for any odd k > 3 provided X k and Y 2k are functionally 
independent. 

This last point can be demonstrated by showing that \X k , y 2k \ 7^ and that such a 
property remains true after projection. As \X k , y 2k \ is a (2k + l)th-order operator, it is 
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enough to prove that the coefficient of one of the highest-order terms, namely D^D^, does 
not vanish. Such a term can be obtained by multiplying the term proportional to D V D 2 in 
one the commutators [Xk,2uHi] by the terms proportional to D 2 in the remaining 2uH,j's 
(j 7^ i) and summing over i, as well as over all permutations p. 
On starting from the relation 
[X k ,2uH Q ] = -[Dl,2uH ] 

rD r 



[r 2 , 2uH }n k + {{rD r ) 2 , 2uH ] + 2 (aK + b) 7Z 2l l, 2ooH 

i 

+ 2(a1l + b)Y^ K 2l Z [rD r , 2cuH ] - cu 2 [r 4 , 2uH ], 



i 

it can be easily shown that the coefficient of D V D 2 in [X k , 2uH ] is equal to 
— 4 sin 99 cos v?- It is then enough to apply TV to this result to prove that the term 
-4 sin (<p + y) cos (ip + y) D^D 2 appears in [X k , 2uHi\. Finally, we get the term AD v D 2k 
with A = —2~ 2h+4 k\ k sin kip cos kip ^ in [^,3^]; thereby showing that Y 2 fc] 7^ 0. 

Taking into account that the truth of the conjecture of p] is already known for k = 1 
13], we conclude that it is valid for any odd k. 

5 The k = 3 case 

To illustrate the approach followed in the present paper, it is worth considering the smallest 
allowed value of k, namely k = 3. In such a case, another proof of the superintegrability 
of Hk is indeed available through the use of other types of differential-difference operators, 
modified boson ones and modified boson oscillator Hamiltonians [13J. Here we plan to 
provide some links between those operators and the corresponding ones constructed in 
sections 2, 3 and 4. 

The Hamiltonian H$ appears when considering the three-body Hamiltonian of the CMW 
problem [6] and eliminating the centre-of-mass motion. The latter Hamiltonian 

#cmw = + u 2 x 2 ) + 2a (-L + ±- + -pj + 6/3 (-L + _L + ^] , 

i=1 \ x 12 x 23 x 31/ \yl2 £/23 #31/ 

where x i; i — 1, 2, 3, denote the particle coordinates, Xjj = Xj — x.,-, z 7^ j, and j/^ = x^ + 
— 2xfc, i j k i, can indeed be separated into a centre-of-mass Hamiltonian H cm = 
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— dx + u 2 X 2 and a relative one H re \, coinciding with H 3 , by defining X = (xi+x 2 + Xs)/y3, 
r cos cp = X12/V2 and r sin cp — y 12 / 

In [13] , Hqmw has been transformed into a .De-extended Hamiltonian Hcmw by realizing 
the 12 operators of Dq in terms of the particle permutation operators Kij and the inversion 
operator I r in relative coordinate space and by introducing three differential-difference 
operators Di, i = 1, 2, 3, defined as in equation (2.5) of [15J. Three sets of modified boson 
creation and annihilation operators a J = (coxi — Di)/\^2co, ai = (coxi + Di)/y/2uj, i = 1, 

2, 3, have then be considered and used to build three modified oscillator Hamiltonians 
hi = {a\, a«}/2, % — 1, 2, 3. In terms of the latter, 

"Hcmw = hi, (5.1) 

i 

while the sixth-order integral of motion of Hqmw has been derived from (2a;) 3 h\ by 
projection in the D$ identity representation space. 

The correspondence <-> TZ 2k+3 l and i^i,. 7?. 2fc Z, where (r/fc) = (123), leads to 
equation (2.6) of [15] . connecting the differential-difference operators of equation (12.21) with 
the ZVs. As a result, we can express the three dependent modified annihilation operators 
Ai, i = 0, 1, 2, defined in (13.81) . in terms of the three independent operators ai, i = 1, 2, 

3, as Aq = (ai — a 2 )/\/2, A\ = — (a 2 — 0,3) /y/2, A 2 = (as — a\)/y/2. The alternative set of 
operators Bi, i = 0, 1, 2, is given by B = (a\ + a 2 — 2a 3 )/v / 6, B\ = — (a 2 + a 3 — 2ai)/\/6, 
-B2 = (03 + ai — 2a 2 )/\^6. From this, it can be easily checked that there is a single relation 
Aq — A\ + A 2 = among the Aj's (and similarly for the B^s) and that all relations in 
(13.13!) and (13. 14p reduce to that one for k = 3. Furthermore, the commutation relations 
(13.151) and (13 . 161) entirely agree with those that can be derived from the above-mentioned 
relations between the A^s and the a^s and from the commutation relations of the latter 
given in equation (4.2) of [IB] . 

Turning now ourselves to the modified boson oscillator Hamiltonians, we get the equa- 
tions 



H = 


\[hi + h 2 - 


\{a\, a 2 } — 


\{a\, ai}], 




H 1 = 


l[h 2 + h 3 - 


\{a\, as} — 


\{a\,a 2 }], 


(5.2) 


H 2 = 


\[hz + h x - 


\{a\, a{} — 
11 


\{a[,as}}, 





where the extra anticommutators on the right-hand sides serve to compensate for the centre- 
of-mass degrees of freedom present in hi, h 2 , h 3 , but absent from H , Hi, H 2 . From 
equations (H~2l) . (IOI) and (|5\2|) . it follows that 



in agreement with the decomposition 

%cmw = H cm + H 3 , H cm = \u{a\ + a\ + a\, a x + a 2 + a 3 } = -d\ + u 2 X 2 . 

Finally, the operator y e , giving rise by projection to the sixth-order integral of motion of 
H 3 , could be expressed in terms of a\, a,i, i — 1,2, 3, through equations (14.41) and ( 15.2ft . 

6 Conclusion 

In the present communication, we have provided an alternative approach to the construction 
of the sixth-order integral of motion of the CMW Hamiltonian, wherein the centre-of-mass 
motion has been eliminated from the very beginning, in contrast with the previous one 
based on the use of differential- difference operators related to the three particle coordinates 
on the line [13] . The present construction instead employs differential-difference operators 
D r , D v on the relative-coordinate space. A key step in the construction has been the 
introduction of three dependent sets of modified boson creation and annihilation operators 
with simple exchange properties under the dihedral group Dq connected with the CMW 
problem. 

An analogous approach, taking advantage of some recent progress made in [15], has 
allowed us to prove the superintegrability of all TTW Hamiltonians H^ with odd k > 3. 
In such a case, the associated group being D 2 k, we have been led to consider k dependent 
sets of modified boson creation and annihilation operators and the k modified oscillator 
Hamiltonians that can be built from them. In terms of the latter, the construction of 
the missing integral of motion Y 2 u has proved easy. Furthermore, Y 2 k has been shown 
to be functionally independent of the first integral of motion X^, thereby completing the 
demonstration. 
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Since for even k values, the operator D v assumes a different form from that considered 
here for odd k, it is not clear yet whether the present formalism can be extended to such 
values and lead to a proof of the remaining part of the conjecture. This is left for a future 
study. 
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